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Abstract 

The study of the propagation of electrons with a varying spinor orientabihty is performed using 
the coordinate transformation method. 

Topological Insulators are characterized by an odd number of changes of the orientabihty in the 
Brillouin zone. For defects the change in orientabihty takes place for closed orbits in real space. 
Both cases are characterized by nontrivial spin connections. 

Using this method , we derive the form of the spin connections for topological defects in three 
dimensional Topological Insulators. 

On the surface of a Topological Insulator, the presence an edge dislocation gives rise to a spin 
connection controlled by torsion. We find that electrons propagate along two dimensional regions 
and confined circular contours. We compute for the edge dislocations the tunneling density of 
states. The edge dislocations violates parity symmetry resulting in a current measured by the 
in-plane component of the spin on the surface. 
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I Introduction 

The propagation of electrons in solids is characterized by the topological properties of 
the the electronic band spinors. Topological Insulators can be identified by 

an odd number of changes of the orientability [8] of the spinors in the Brillouin zone. As a 
results non trivial spin connections with a non- zero curvature characterized by the Chern 
numbers can be identified. In time reversal invariant systems one finds that for Kramer's 
states the time reversal operator T obeys = — 1 and one thus the second Chern number 
for four dimensional space is given by (—1)'^ = —1, where u is an odd number of orientability 



changes 



24| 



Real materials are imperfect and contain topological defects such as dis 



18| ,disclinations [l9|, l20| and gauge fields induced by strain in graphene |2l|, |22| ;therefore, a 



ocations 



15 



natural question is to formulate the physics of Topological Insulators in the presence of such 
defects |8|. These topological defects can be analyzed using the coordinate transformation 
method given in ref.|26| which modifies the Hamiltonian for a Topological Insulator with a 
defect by the metric tensor and the spin connection 30l-l34l|. 

The effect of strain fields dislocations and disclinations plays an important role in material 
science and can be study using Scanning Tunneling Microscopy {STM) and Transmission 
Electron Spectroscopy {TEM ). Therefore we expect that the chiral metallic boundary 29 1 
will be sensitive to such defects. 

In t his p aper we will introduce the tangent space approach used in differential geometry 



24 



33 



34| to study propagation of electrons for a space dependent coordinate |26[. We find 



;hat the continuum representation of the edge dislocation 26| generates a spin connection 



30 



32| which is controlled by the Burger vector. 



Using this formulation we obtain the form of the topological insulator in three dimensions 
which simplifies for the surface Hamiltonian (on the boundary). For the surface Hamiltonian 
we find that the electronic excitations are confined to a two-dimensional region and to a set 
of circular contours of radius Rg{n). 

The contents of this paper is as follows: In chapter // we introduce the gemetrical method. 
In section II A we present the geometrical method for the edge dislocations and strain 
fields. In section II B we consider the effects of the strain fields on the three- dimensional 
Topological Insulator (TJ). The Chiral model for the boundary surface is presented in 
section III A. Section IIIB is devoted to the derivation of the metric tensor and spin 
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connection for an edge dislocation [26|]. In section II IC we identify the stable solutions. 
Section HID is devoted to the stable two dimensional solutions n = and section HIE is 
devoted to the stable solution for circular contours n = ±1. Chapter IV A is devoted to the 
computation of the tunneling density of states. In section IVB we present results for the 
two dimensional region n = 0. Section IVC is devoted to a large number of dislocations. In 
section IV D we compute the tunneling density of states for the circular contours n = ±1. 
In chapter V we consider the current which is given by the in-plane spin component. In 
section VA we show that this current is zero for a TI. In section VB we show that in the 
presence of an edge dislocation the parity symmetry is violated, and current, representing 
the in- plane spin component, is generated. Chapter VI is devoted to conclusions. 

II-The Geometrical method for dislocations and strain fields 

A-General Considerations 

A perfect crystal is described by the lattice coordinates f = [x,y,z]. For a crystal with 
a deformation , the coordinates r are replaced hjf-^R = f+u = [X^(r), X^(r), X^(f)] 
where u{r) is the local lattice deformation and X"-, a = 1, 2, 3 is the local coordinate which 
changes when we move from one point to another. 

In a deformed crystal we introduced a set of local vectors which are orthogonal to each 
other (cb, Ca) =< e^|ea >= 5^ and local coordinates a = 1, 2, 3. The unit vector Ca can be 
represented in terms of a Cartesian fixed frame space with the coordinate basis ,/i = 
x,y,z. In the fixed Cartesian frame the coordinates are given by x'^. Using the Cartesian 
basis we expand the deformed medium in terms of the local tangent vector Ca : Ca = e^9^ 
(for the particular case where vectors Cq are given by Ca = da, the transformation between 
the two basis is = 5^). Any vector X (in the deformed space) can be represented in 
terms of the unit vectors Cq or the (9^ (the tangent vectors in the Cartesian fixed coordinates 
space). The vector X can be represented in two different ways, X = A^e^ = X^d^ (when an 
index appear twice is understood as a summation, A^Ca = X]a=i 2 3 "^'^^a)- The dual vector 

is a one form and can be expanded in terms of the one forms dx^. We have: = e^dx, 
where represents the matrix transformation e° = {df^X"')dx^. The scalar product of the 
components e^e" = g^^^, ej^e^ = 6a^b defines the metric tensors, g^^^, (in the Cartesian frame 
) and 6a b in the local medium frame. 
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B- Application to the Topological insulators in three dimensions 

The three dimensional electronic TI bands for Bi2Se3 and Bi2Te^ can be represented 



using four projected states [l6j, \orhital = 1, 2 > ®\spin =t, \.> (the Pauli matrix r describes 
the orbital states and the Pauli matrix a describes the spin). The effective h^^ Hamiltonian 
in the first quantized form is given by: 

h^^ = hvo[ky{ai (g) ri) - h{(j2 ® n) + ek^{a^ ® n) + M{k){I ® rg)] (1) 

The parameter M{k) determines if the insulator is trivial or topological. For Bi2Se^ and 
SisTes the gap is inverted^amelvM {k) = -MQ+Bikl+B2k\ with Mq > 0, 5i > 0, ^2 > 
and therefore topological [111, Q, 16 1. 



Using the metric tensor g^^y given by the coordinate transformation ( the transformation 
between the two sets of coordinates - the one without the dislocation and the second with 
the dislocation ) e^e" = g^^y, defines the Jacobian a/G where G = det[g^^y\. We find that 
the derivative for a spinor component \E'(°)(r), a = [1 = 1 t; 2 = 1 J,; 3 = 2 f; 4 = 2 |] is 



replaced by the covariant derivative 



3l|: 



V^vl/(")(f) = d^¥''\f) + ^J;'''\T'', t%¥^\r} (2) 
where f'' ,a = 1,2,3,4,5 are the matrixes: = -P^ = -(as ® n); = pi = (ai (g) ri); 

f3 = r3 = ((T3 ® Ti); = P4 = (/ ® T2y,P = P^ = (/ ® Tg). 

The spin connection cu^'^ determines the covariant derivative 
tangent vectors e°: = d^X°-{f); a = 1, 2, 3 ; fi = x,y,z. 



3l\ is given in terms of the 



_leP."e-'^(9,e.,,-9<,e,,,)e^ (3) 

We notice the asymmetry between e'"'" and Cay. e^''^ = g^'^e^ and ea,u = ^a.feC^. As a result 
the Hamiltonian in eq.(l) in the second quantized form is replaced by: 

i/O^) = hvoj d^rVG[^^{f)[e^t''{-iV^) - Ef{I ® I) + t%-Mo)]'^{r) 

+5i(7'^''^(V^^I(f))(V.^i(f)) - Big^^'{V^^l{^V,^2)] 

(4) 

where e^it" = J^a^^a'^" = ^^i^), [f^(f)f^(f) + f^(f)f'^(f)] = 2^'^''^(f) , rfet[^^'^(f)] = G 
and is the covariant derivative given in terms of the spin connection given in equation 
(2): V^^(")(f) = 9^^(")(f) + lJ;''\t^,t%^^^\r) 



C-The Mechanical strain effect on if^^D) 



From the work of p/7[ we learn that the effect of the strained field is different on Bi2Se3 
than on Bi2Te^. In Bi2Se?j the compressive strain decreases the Coulombic gap while 
increasing the inverted gap strength induced by the spin-orbit interaction. We will use the 
result in equation (4) to analyze the effect of strain. The strain field ejj (symmetric in 
is related to the stress field cxij and elastic stiffness Lame constant A and /x: cTij = 
^^i,j^k,k + 2/uejj-. Applying a constant stress CTjj one can determine the value of the constant 
strain field tij which is related to the tangent vectors e* = 5ij + e^j-. In the present case the 
spin connection and the Christofel tensor vanish. The metric tensor gi j is given by :gi^j = 
6ij + 2eij. Using this formulation we can investigate the effect of the stress on the Bi2Se-j, at 
the r point = 0. The Tl Hamiltonian given in eq.(4) ,M{k) = -Mq + Bikl + ^2(^2 + kl) 

n 

with the inverted case Mq > [12] . The Hamiltonian in eq. (4) is replaced by: 

H^^^-^'"''^'^^ = hvo J rfVv^[M/t(f)[f'^(5^,, + e^,J(-29^)) + f^(-Mo) + 5(1 - 2e^,,)9/U]v[/(f) 

^ hvo J d\VG[^^ir)n5^,a{l+ < e >)i-td^))] + f^(-Mo) + 5(1 - 2 < e > 5^,,)d^t^d,]^{r) 

(5) 

In equation (5) we have used the average strain field < e >, < e >= '^i-^+'^^.^+'^a.s ^ replace 



the spinor field '^{f) by ^(^)a/(1+ < e >) = ^(r). As a result we obtain: 

^(3D-strain) ^ ^ f d^rVCl^^f^T^' {-id.) + T \ ^"^"^ , + 5 ^ ^ " ^ ^ ' ^ fU] ^(r ) 

(6) 

For the compressive case < e > is negative, < e >= — | < e > | . As a result we observe 
that the inverted gap is enhanced = (ii|'<°>|) > l^ol- 

In the same way we can show that the Coulomb interaction is reduced: We introduce the 
Hubbard Stratonovici field Oq to describe the Coulomb interactions. 

H'-' = J rfVv^[/(-e)-aoVl/t(f)vI/(f) + il^|^ao9^9,ao] 
^ J £rVG[I{-e) ■ aoVl>t(f)vl/(f) + ~ ^ ' ^K od^d^ao] 



(7) 



Next we rescale an = , and obtain: 



J V 1 — 2 < e > 
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We observe that for the compressive case the effective charge Ceff. = ^i^2<e> ~ y/i+2\^ > | 
reduced and therefore the Coulomb gap decreases, while at the same time the inverted gap 
increases, = (ii|^<°j>|) > l^ol in qualitative agreement with 17 1. 



III-The chiral metal with an edge dislocation 
A-Description of the Chiral model 

The low energy Hamiltonian for the bulk 3D TI in the Bi2Se^ family was shown to 
behave on the boundary surface (the x, y- plane) as a two dimensional chiral metal j?! . 

(9) 

h^'^ = hvF[i<y^dy — ia'^dx] is the chiral Dirac Hamiltonian in the first quantized language. 
vp ~ 5-10^^ is the Fermi velocity, cr is the Pauli matrix describing the electron spin 
and /i is the chemical potential measured relative to the Dirac F point. The Hamiltonian 
for the two dimensional surface L x L describes well the excitations smaller than the bulk 
gap of the 3D TI at 0.3 eV. Moving away from the F point, the Fermi velocity becomes 
momentum dependent; therefore, we will introduce a momentum cut off A to restrict the 
validity of the Dirac model. The chiral Dirac model in the Bloch representation takes the 
form: h = hvpiK x a) ■ z = hvF{—cr^ky + cr'^k^) The eigen-spinors for this Hamiltonian are 
: \u{K) >= [\u^iK) >, \ui{K) >f = \K > ze*>^('=-'=^)]^ where x{kx,ky) = tan-\!^) is 
the spinor phase and e = hvp^yk'^ + ky is the eigenvalue for particles . For holes we have the 
eigenvalue e = —hvp^Jkl + fc^ and eigenvectors \v{K) >= [Iv-^i^K) >, \vi{K) >]'^ = \K > 
ie^^^'^'^'^'y^]'^ . The chirality operator is defined in terms of the chiral phase x{kx, ky): 

K 

(<? X — ^) ■ z = sm[x{kx, ky)]a^ - cos[x(A;^, ky)]a^ (10) 
\K\ 

The chirality operator takes the eigenvalue — (counter-clockwise) for particles 

[sm{x{kx,ky))a^ - cos{x{kx,ky))a'^]\K > ®[1, ^e^^(^-^'')]^ = -\K > ze^^C^-'^^)]^ and 
+ (clockwise) for holes [sin(x(/cz, fc?/))cr^ - cos{x{kx, ky))a'^]\K > (g)[-l, ■ie*^^''^'''^)]^ = \K > 

B-The effect of edge dislocation on a two dimensional chiral surface Hamilto- 
nian 
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We use the notation x'^ ,fi = x,y and X" ,a = 1, 2 to describe the media with dislocations. 
For an edge dislocation in the x direction the Burger vector B^'^^ is in the y direction . The 
value of the burger vector B^"^^ is given by the shortest translation lattice vector in the y 
direction. (For the TI Bi2Se-i the length of the vector B'^'^^ is 5 times the inter atomic 
distance ). Following [26.] we introduce the coordinate transformation for an edge 
dislocation: r = {x,y) — )■ [X{r) = x,Y{r) = y + tan^-'^(|)] with the core of the 
dislocation centered at r = (0,0). The matrix elements fields for the edge dislocation is 
given by : 

e^ = 9^X"(r); a = 1,2; fi = x,y (11) 
We express the Burger vector in terms of the the partial derivatives with respect the coor- 



dinates a = 1, 2 in the dislocation frame and fi = x,y for the fixed Cartesian frame 2q |: 



S.4 - = B'"S'( f) (12) 

Using Stokes theorem, we replace the Hne integral /d.'-ej(r) by the snrfaoe integral 

J J dx^dx'^[dxey — dye1]. For a system with zero curvature and non zero torsion T^fJ we find 
that the surface torsion tensor integral J J dx^dx^Tj^)^ is equal to j J dx^dx'^ld^e^ — dye"^], 
and therefore both integrals are equal to the Burger vector. 

(13) 

where dx^dx'^ represents the surface element. The tangent components can be expressed 
in terms of the Burger vector density B^'^^5'^{r) 



m ■ 



2n ' (a;2 + y2) ' v ^ ^-k ' (x^ + y'^) 

el = l- ej = (14) 

Using the tangent components, we obtain the metric tensor g^^y. 

elel ^ e\el + eje^, = g,A^- e^e^ ^ ele\ + e^ej = 5,,, (15) 

The inverse of the metric tensor gp^^vir) is the tensor g^'^ir) defined trough the equation 
9pu,T{f)9^'"{f) = Using the tangent vectors, we find to first order in the Burger vector 



the metric tensor (7^ jy and the Jacobian transformation v G: 

_B^'^ y B^) y _ 5(2) y 

9x,x J-) 9x,y „ 2 I 2' ^y^y n 2 I 2' ^l''^ ' |.fl'A',i^J -I- 



(16) 

The inverse tensor is given hy.g^'^ ~ 1, g^'^ = g^'^ = 2I. 2 ; S'^'^ = 1 + tt-j- Using 

the inverse tensor g^^^^ we obtain the inverse matrix which is given by: 

= ea^ug"'^ = {6a,tel)g'''^ = etg'''^ (17) 

Using the components we compute the the transformed Pauh matrices. The Hamilto- 
nian in the absence of the edge dislocation is given by h^'^' = i'~i°'da = J2a=i 2 ^7"^a where 
the Pauh matrices are given by 7^ = — cr^ ^^2 = 0"! and 7^ = a^. (We will use the conven- 
tion that when an index appears twice we perform a summation over this index.) In the 
presence of the edge dislocation, the term 'j"'da must be expressed in terms of the Carte- 
sian fixed coordinates fi = x,y. As a result, the spinor \E'(r*) transforms accordingly to the 
SU{2) transformation . If \E'(-R) is the spinor for the deformed lattice, it can be related 
with the help of an SU{2) transformation to the spinor \E'(r*) in the undeformed lattice: 
\E'(X, y) = e~* ^V^'°"^vl>(x, y) . Where Sip{x,y) is the rotation angle between the two set of 
coordinates: 5ip{x,y) = tan~^{Y) ~ tan~^(^). Using the relation between the coordinates 
X = X, and Y = y+ ^^tan^^{^) with the singularity at x = y = gives us that the deriva- 
tive of the phase which is a delta function, dxSip{x,y) = —dy6ip{x,y) oc 6'^{x,y). Combining 
the transformation of the derivative with the 5*0(2) rotation in the plane, we obtain the 
form of the chiral Dirac equation in the Cartesian space (see Appendix A) given in terms of 



the spin connection 10]^ 24| : 



Z7a9.vl/(i?) = tba^.^^'d^R) = tre^^id, + ^[7^ 7>S*(r1 (18) 

The Hamiltonian h'^-^- — )■ /i^'^f^ is transformed to the dislocation edge Hamiltonian with the 
explicit form given by: 



(19) 



To first order in the Burger vector we find : u]^ = —oj"^ = and —UJy^ = cu^'^ = — -^^^^(r), 



see eqs. (72 — 74) in Appendix A. 



(20) 



In the second quantized form the chiral Dirac Hamiltonian in the presence of an edge dislo- 
cations is given by : 



(21) 



fledge jg ^.j^g Hamiltonian in the first quantized language, /i is the chemical potential and 
\l'(r) = [\E'-f-(r*), \['j^(r)]^ is the two component spinor field. 

C- The Identification of the physical contours for the edge Hamiltonian h'^'^^^ 

In order to identify the solutions, we will use the complex representation. The coordinates 
in the complex representation are given hj, z = ^{x + iy) , 1 = \{x — iy), dz = dx — idy, 
dz = dx + idy. In this re pre sentation the two dimensional delta function 5'^{r) is given by 
5\r) ^ 19.(1) = 



36 



37I . We will use the edge Hamiltonian h^'^^^ and will compute 
the eigenf unctions u^{z,^) = [U^^{z,z),U^i{z,z)]'^ and v^{z,z) = Vt,^{z,z),V^i{z,z)]'^. The 



eigenvalue equation is given by: 



eU,^{z,z) = -[dz + {^)d.{h]Uei{z,z) 
y Ztt z 



(22) 



u^{z, z) = M{z, z)F^{z, z) = 



e 27r 



The eigenfunctions u^{z,z) and v^{z,z) can be written with the help of a singular matrix 
M{z,z) y : 

' Ol fF^,iz,z) 

Fei{z,z) 

(Fe(2;, ^) and F_^{z.,^) are the transformed eigenfunctions for e > and e < respectively 
.) In terms of the transformed spinors the eigenvalue equation h'^^^'^{z,z)u^{z,z) = m<^{z,z) 
and F^\^{z,'z) becomes: 



13(2) 

e"— 



( F^{z,z)\ 




\F.i{z,z) ) 





I{z,z) 




-dz 


( F^^izz) \ 


{i{z,zy _ 




dz 


\Fei,iz,^) I 



_ b(^) /z-z\ n B(^) , iy _ q bC^) / -iy _ 

where 7(2;, z) = e"^^*-^^ = e ^^+^ , (/(z, z))* = e , 1 1(2;, z)| = 1. We search 

for zero modes e = and find : 

d,F,^iz,z) = d^F^iz,z) = (23) 

The solutions are given by the holomorphic representation F^=o^{z,z) = and the anti- 
holomorphic function ^^=04,(2^)^) = fii'^)- The zero mode eigenfunctions are given by : 

u,=oM = e-^(^)/^(z), u,=o,i(z) = e-^(i)/^(^) (24) 

Due to the presence of the essential singularity at 2; = it is not possible to find 
analytic functions /^(-z) and fii'z) which vanish fast enough around z = such that 
J d"^ z{ue=o,\{z))*Ue=o,x{z) < oo. Therefore, we conclude that zero mode solution does 
not exists. The only way to remedy the problem is to allow for states with finite energy. 
In the next step we look for finite energy states. We perform a coordinate transformation 

z^W[z,z]; z-^W[z,z] (25) 

We demand that the transformation is conformal and preserve the orientation. This restricts 



the transformations to holomorphic and anti- holomorphic functions 36|]. This means that we 
have the conditions dzW[z,z] = and dzW[z,z] = 0. As a result we obtain W[z,z] = W[z] 
and W[z,z] = W[z], which obey the eigenvalue equations: 

eF,^{W,W) = -dwF,^{W,W) 
eF,^{W,W) = d^F,^{W,W) 

(26) 

This implies the conditions = {I{z,^))* and ^^^^ = I{z,^). Since I{z,^) is neither 

holomorphic or anti-holomorphic and satisfy |J(2;,z)| = 1, the only solutions for W[z] and 
W[z] must obey I{z,z) = 1: 

J(z,^) = e^^^^^ = e^^™; n = 0, ±1, ±2.... (27) 

For I{z,z) 7^ 1 one obtains solutions which are unstable . The stable solutions will be given 
by a one parameter s curve (s is the length of the curve) r(s) = [x{s),y{s)] which obey the 
equation I{z,z) = 1. The curve r{s) allows us to define the tangent t{s) and the normal 
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vectors N{s). This allows us to introduce a two- dimensional region in the vicinity of the 

— * — * 

contour of r(s) — >■ R{s, u) — f{s) + uN{s). 

IID- The wave function for the edge dislocation-the n = contour 



The condition I{z, z) = e ^ ^x^+y'^' = 1 for n = is satisfied for y = and large value of 



y which obey 2-^(^5^) << 1 ■ The values of y which satisfy this conditions are restricted 
to I{z, z) — e ^x^+y^' Ki 1. This condition is satisfied for values of y in the range: 



2^(Z^)<^<T<1 (28) 



y TT 

We introduce the radius Rg — ^ and find that the condition I{z,^) 1 gives rise to 
the equation for y. The solution is given by + (y ± ^Rg)^ — (^RgY- Therefore, for 
\y\ > \d\ > {^)2Rg > 2Rg we have I Ki 1 which corresponds to a free particle eigenvalue 
equations. 

eFet(a;, y) = e ^ (-'+y'') [-d^ + idy]F^i{x, y) 
a; [-d^ + idy]F^^{x,y)- 
eF^{x, y) = e - % + idy\F,^{x, y) 

^ [d^ + idy]F,^{x,y) 

(29) 

For \y\ > d the eigenfunctions are given by: Ue,'[{x,y) = e (x, y), Ue^i{x,y) = 

_s(£)./^_^ 

e ^x-iy) p^^^(^x,y) where F^'^{x,y) and F^{x,y) are the eigenfunctions of equation (21). 
The envelope functions e ^x+ij/''^ g ^tt ^^-ij,'' which multiply the wave functions F^^[x, y) , 
Felix, y) impose vanishing boundary conditions for the eigenfunctions Ue,i{x, y) and Ue,-^{x, y) 
at y — >■ ±oo . therefore, we demand that the eigenfunctions Ue^^^ix, y), U^^iix, y) should vanish 
at the boundaries y — Since the multiplicative envelope functions for opposite spins is 
complex conjugate to each other wc have to make the choice that one of the spin components 
vanishes at one side and the other component at the opposite side. Two possible choices 
can be made: 

b(2) , 1 ^ _B^, 1 . 

C/,,^(x, y = f ) = e '^'Fe^ix, f ) = C/,,^(x, y = -f ) = e ^^^^ V^(x, -f ) = 
or 

b(2) . 1 . b(2) , 1 s 

C/,,t(x,y = -f ) = e ^^+^(-*) V,t(x, -f ) = C/a(x,y = f ) = e ^'^ F^(x, f ) = 
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Making the first choice, (both choices give the same eigenvalues and eigenfunction) we 
compute the eigenfunctions Ff^^{x, y) and F^{x, y) for \y\ > d using the boundary conditions 

Fetix, y=^)=Q; F^{x, y = ~)=Q (30) 

Due to the fact that the solutions are restricted to \y\ > d no conditions need to be im- 
posed Sii X — y — In the present case we consider a situation with a single disloca- 
tion. This is justified for a dilute concentration of dislocations typically separated by a 
distance / lO^^m. ( In principle we need at least two dislocations in order to satisfy 
the condition that the sum of the Burger vectors is zero.) The eigenvalues arc given by 



e = ±hvF-\Jp^ + q^. The value of p is determined by the periodic boundary condition in 
the X direction p(m) = = m = 0, 1, {N — 2), [N — 1) and a is the lattice 

constant a ~ The value of q will be obtained from the vanishing boundary conditions at 
y = ±1^. The eigenfunctions F^^^{x,y) will be obtained using the linear combination of the 

spinors introduced in chapter ///. In the Cartesian representation we can build four spinors 
rp^q{x,y), rp^_q{x,y),r_p^q{x,y),T_p^_q{x,y) which are eigenstates of the chirality operator 
and are given by: 



1 




r,,,(a;,|/)=e^^V«^ 



rp^^q{x,y)^e'^^e 



_,,_,(x,y) = e-^^e-^^| \ \ (31) 

where tan[x{p,q)] = |. 

The TI Hamiltonian h^-^{x,y) — hvp[i(T^dy — ia^d^] is invariant under the symmetry 
the operation x — > —x which is described by the transformation : PxxP~^ — —x] 
P^a'P-' = -a'; P^yP'^ = y; P^a^P'^ = a\ 

The edge Hamiltonian ff^^^ contains in addition the term a'^5{r) which changes sign under 
the symmetry operation P^ . As a result the symmetry operation does not commute with the 
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edge Hamiltonian [/i^'^^'^, P^] 7^ 0. This result demands that we construct two independent 
eigenfunctions F^^^^^\x,y) {right — mover) for p > and F^p^Q^\x,y) {left — mover) 
p<0. 

FiloT\^: y) = ^(?)rp,,(x, y) + B{q)T^,_,{x, y) 
F-pS;!^\x,y) = C{q)T_p^g{x,y) + D{q)T_p^_g{x,y) 

(32) 

Employing the boundary conditions given in equation (29) we obtain the amplitudes ^||y , 
and the discrete momenta q+. Using the pair rp^q{x,y) , Tp_q{x,y) p > we obtain : 

e(p>0,9+),4 



'x, y) = ie^f^e5>^(*''«+)[e^(«+^+^^(^''«+)) + (-l)^+ie-^(«+2'+5x(p,<z+))]. \y\ > ^ 



g = g+ = y/c + y tan ^(— ); k = l,2, 3...; tan[x(p, q+)] = ( — 
L L p p 



271 

^ip, Q+) = ±hvFd (-^m)2 + ql 

(33) 

Similarly, for the second pair F _p^q{x , y) ,T _p^_g{x , y) , p > we obtain: 

Fei-p>o\.)i(^^y) = -ie-^P^e-5>^(*'''?-)[e^(«-^-^x(p,?-)) + (_i)fe+ie-i(9-y-|x(j',?-))]. |^| > 

q = q^^yk- y tan-^(— );A; = 1, 2, 3...; tan[x(p, ?-)] = (— ) 
L L p p 



e{-p, g_) = ±hvF\j {^m)'^ + ql 



For the state with zero momentum p = we find: 



(34) 



^r.:o2.t(^> y) - 2e-t^ cos[gy +^];\y\>d 



^rP=M,,(^, y) - ^2e^ cos[gy - ^]; |y| > d 
? = ^ + ^A;;A; = 0,l,2,3... 
e{p — 0,q) — ±hvF\q\ 



(35) 
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The eigenfunctions for the dislocation problem for \y\ > d will be given in terms of the 
envelope functions e ^tt '-x+iyi ^ g ^vr '^x-iy' (U^.^(^x,y) — e ^vr ^'^+iy' Ff^^[x,y), U^^_i_{x,y) — 
e 2- ^^-iy'F,^i{x,y)). 

The explicit solutions are given by : 

u't=''''\x,y) ^ [U^;;=''\x,y),uil=''''\x,yr; ut='''\x,y) 

[U^^=''\x,y),ut'''\x,y)f. 

The components of the spinor are given by: 

t^t y) ~ Gi(a;,|/)L ' ^e(-p>o„_),t(^' y) 

^ ^ '^^ GHx,y)L ^ei-p>o,,.u^^^y) 



(36) 

where G(x,y) = 1 — ^^|^ y2(J+;/2) Jacobian introduced by the edge dislocation. The 

eigenstates are normalized and obey: J dx J dy^jG{x, y){U^~^"^\x, y))*U^~^'^\x, y) 
Sa,a', J dx J dy\/G{x,y){u'^~^'^\x,y))*U^J^~^'^\x,y) ^ 5o-,o-'. The normalization factor 
2const^B( )) ^ -^^^ ^ wcak dcpcndcncc on the Burger vector B^'^^ . This dependence is a 
consequence of the Jacobian ^/G which affects the normalization constant. 

(The multiplicative factor e ^tt yx±iy> gives rise to a weak non-orthogonality between the 
states. This non-orthogonality of the linear independent eigenfunctions can be corrected 
with the help of the Grahm-Shmidt method.) 

For the present case, backscattering is allowed but it is much weaker in comparison to 
regular metals. This is seen as follows: Time reversal is not violated; due to the parity viola- 
tion, the eigenstates ui"'~^'^\x, y) ,Me'^~°'^^(a;, y) are not related by a time reversal symmetry 
{Tu^r~^'^\x^ y) 7^ u^r~^'^\x^ y)) . As a result, the backscattering potential T^j.-p controlled 
by a finite matrix element between states with different eigenvalues e(— p, q^) ^ e{p, g+) (con- 
trary to regular metals where the impurity potential 14>,-p connects states with the same en- 
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ergy). In the present case |e(— p, q-) — q+)\ — hvpHy (x"^)^ + Q-" y (x"^)^ + ^+11 ^ 
the eigenvalues are not equal, therefore the finite matrix element controlled by the backscat- 
tering potential gives rise only to a second order backscattering efi^ect! 

HIE- The circular contours-the wave function for n 7^ 

2 ^^^^ ( '^t") ) 

The equation I{z,^) = e ^^ITjT^^ = gj27m gjygg ^j-^g Qf ^.j^g contours for n = 



±1,±2,±3, ... shown in figure 1. The radius Rg for the fundamental contour(n = 1) is 

Rg 

H' 

{x{s)r + {y{s)±Rg{n)r^{Rg{n)r (37) 



represented in terms of the Burger vector B^'^\ Rg — and Rg{n) — 



The centers of the contours are given by :[x^y\ = [Q,Rg{n)] for n ^ 0. When n > the 
center of the contours has positive coordinates (upper contour) and for n < the center 
has negative coordinates (lower contour). Each contour is characterized by a circle with a 



radius Rg{n) = ^ centered at [x — 0,y — Rg{n)]. The contour is parametrized in terms 
of the arc length < s < 27r-j^ which is equivalent to < < 27r . Each contour is 
parametrized by r(s) = [x{s), Rg{n) + y{s)] where x{s) = -Rg(n) cos [-^^^] = Rg{n)cos[Lp] 
and y{s) = Rg{n) sin[-^^^] = Rg{n) sm.[ip\. We will extend this curve to a two dimensional 
strip with the coordinate u in the normal direction: For the curve curve r(s) — [x[s),y[s)] we 
will use the tangent t{s) and the normal vector N{s) Therefore, the two dimensional region 
in the vicinity of the one parameter curve r{s) is replaced by r{s) — >■ R{s, u) — f{s)+uN{s). 

x{s,u) — Rg{n) cos[—^—-] + ucos[- 



Rg{n) Rgin) 
,(.,«) = i?,(n)sin[^]+«sin[^] 

(38) 

We will restrict the width \u\ such that e*27rng±j?j j-^. where r] obeys < | < 1 , < 
— Rg{n) ft! Rg{n){-^) < "^g^"-* . In these new coordinates, the Dirac equation is 

27rn 



approximated for \u\ < Rg{n){-^) — by : 

eF^{s,u) = -I{s,u)e-'^[d^ - \ ds]F,^{s,u) ft -e-*W[a„ - ids]F,^{s,u) 
eF,^{s,u) = {I{s,u)ye'^[d^+ ^ \ ds]F,^{s,u) ft e^[d^ + ids]F,^{s,u) 

^ + Rain) 

(39) 
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The solution for the contour n 7^ 0, < s < 2nRg{n); \u\ < 

The periodicity in s allows us to represent the eigenfunctions in the form: F^^{s,u) = 
Er=-ooE,e^'^*^e'^-F,t(j,g) and F^{s,u) = E,°l-oo e^^^'^'^^^^e^^^F^^ g). We 
find: 

Kg[n) 



j + 1 
Kg[n) 



(40) 



The determinant of the two equations determines the relation between the eigenvalue e, the 
transverse momentum Q{e) and the eigenfunctions F^{j,q),Ff^{j,q). The eigenvalues are 
degenerate and obey : e(j — l;k) — e(j — —{I + l);k) ,where I > 0. 



±g(e); Q{e) = Je 



2 _ f ^ ~^ 2 \2 



F,(l, q) ^ [F,^(l, q), F,^(l, q)f oc [1, g-^'^W)]^; k{Q, I) = tan'^^^ 



(41) 



The value of the transversal momentum Q{e) will be determined from the boundary con- 
ditions at We will introduce a polar angle 9 measured with respect the Cartesian 
axes: The angle < ip{n = 1) < 27r for the upper contour n = 1 centered at \x = 0,y = Rg] 
is described by the polar coordinate < ^ < tt measured from the center of the Cartesian 
coordinate [0, 0]. The lower contour centered at [x — 0,y — —Rg] characterized by the angle 
< (p{n — —1) < 27r is described by the polar angle 9 restricted to tt < ^ < 27r. We 
establish the correspondence between (fi{n — ±1) and 9: 

(fi{n — 1) —29 + — for the upper contour n — 1, < 9 <7r 

(^(n = —1) = 26* + — + TT jar the lower contour n = —1, < ^ < tt 

(42) 

Following the discussion from the previous chapter we will introduce the following boundary 
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conditions: 

Ft"\^^^ = ) = 0; F^r'\'^y = = 

D{n ^±1) = D (43) 

For the two contours n = ±1 we introduce eight spinors T[^'f^'^^\ip{n = 
±l),u),rl'rQ'\^{n = ±1),k), rLtc?'n<^(^ = ±1),«), rLt_Q^((^(n = ±l),u). Using this 
spinor we will compute the eigenf unctions. For the case n = we had only four spinors 
given in equation (30). The four spinors have been used to construct the eigenfunctions 
Fp>o^q^^ (x, y) for p > and Flp~Q (x, y) . Due to the fact that for each n 7^ we have two 
contours n = ± we have eight spinors which will be used to construct the eigenfunctions. 



rlTo'\<P{n - ±1),«) = e*'(^("=±i))e-*«" 



1 

gi(v3(n=±l)) g-jK(i,Q) 
1 

^i{<f{n=±l)) ^iK{l,Q) 
1 

_g-»(v("=±i))g««('>Q) 



r^T o'H<^(^ = ±1), = e-^'(^("=±^))e-^«" I ^ 1 (44) 

-l,-Q VrV / I _g_i(^(„=±i))g_i„(,,Q) I ^ ^ 

Using the vanishing boundary condition given in equation (42) we construct for this case 
similar spinors as the one given in equation (31). In the present case we have for each 
n ^ two contours, therefore the number of spinors will be doubled. We find instead of the 
eigenfunction given in equation (33) two sets of eigenfunctions with momentum Q_ (which 
replaces , see (33)) and (which replaces g+ , see (32)) . 

Using the boundary conditions given in eq.(35) we determine the quantization conditions 
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Q-,Q+ and the eigenf unctions for the n — 1 and n — —1 contours. 

Q- = ^A; - i tan-^(^^^), k^l,2, 3...; tan[/.(i, Q_)] = 

Q^-^k + ^ tan-^(^±^), = 1, 2, 3... 
tanK/,g+)] = (^±^) 



e(Z,g+)=±n^;^^/(-^)2 + Q^ (45) 

Using the fact that the combined wave function on the contours n — 1 and n — ±1 must 
be finite we obtain two sets of wave functions. We include the envelope function and obtain 
the wave function for Q_ and Q^: The envelope functions e ^tt ^x+iy) ^ g ^tt yx-iy> when 
projected to the contours take a complicated form. The envelope functions can be expressed 
in terms of the functions r]{u) and ({9,u): 

_ Rg{l) \U\ 

_5(2) 



2TT{Rg{l) + 'u)((sin[2^])2 + {rj{u) - cos[2^])2) 

(46) 
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We find for g_: 

+ (_l)'e-^('^'")^"^P^]e-^^(^'«)(-^'(«)+^°^l2ei)e-*^(2^+^) sm[Q_u - -k{1, g_)]]; 



(47) 



Similarly for g+we obtain the wave function: 
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^gC(e,n)sin[2e]g-^C((^,n)(^(n)-cos[2e])g-^^(2e+f ) cos[g+K - ^/^(i, g+)]]; 

U^a,Q^liO:u) = 2ze^'^(''«+)[-(-l)^e-^('^'")^^^[2e]g.c(M(-r,(n)+cos[2e])g^(/+i)(2^?+f )3i^[g^^ _ g^)] 

^gC{e,«)sin[2e]giC(e,«)W«)-cos[2e])g-i(«+l)(2e+^) gj^Jg^^ ^ g^)]]. 

UeiiTui^^ ^) = 2e^''(^'«+)[-(-l)^e-^(^''^)^'"[2e]g<(^?,«)(-,W+cos[2e])g^a+l)(2^?+f ) ^^^^g^^ _ 1^^^^ g^)] 

^gC{e,«)sin[2e]giC(e,«)W«)-cos[2e])g-i(i+i)(2e+^) cos[g+« + Q+)]]; 



1 

2' 

(48) 



where G ^ {9, u) is the Jacobian transformation induced by the metric tensor. 
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IV -Computation of the STM density of states 



A-Description of the STM method 

The STM tunneUng current / is a function of the bias voltage V which gives spatial 
and spectroscopic information about the electronic surface states. At zero temperature, 
the derivative of the current with respect the bias voltage V is given in term of the single 
particles eigenvalues: e(m, g_) = ±hvF\/ (^m)^ + g^, e(m, g+) = ±.hvF\ {^rnY + ,111 = 



0, 1,2,3... for contour n = 0. For the upper and lower circular contours n = ±1, we have 
:e(/,Q_) = ±hvF^i^y + Ql ,e(/,g+) = ±hvF^ + 0% ,1 = 0,1,2,3... The 
STM density of states is computed for a voltage V between the STM tip and the sample. 
The tunneling current is a function of the bias voltage V and the chemical potential fi > 



23|: 



oc D{E = eV- s,u) = Y^ D^'^^E = eV; s, u) 



dV 



EE E Ei^'^"^°^"'''^^(^'2/)i'^[e^+/^-^^^i^\/(^^)' + 9.'] 



ri=± m qr=q+,q~ 



EE E Eif^^"=^''''''H^,«)p5[e^+/^-^^^Ft/^^)'+^^ 



n=±l 



(49) 



[f] = + corresponds to electrons with energy < e < /i and rj = — corresponds to electrons 
below the Dirac point e < 0. For the rest of this paper we will take the chemical potentials 
to be /i = 120m (this is typical value for the TI ). We will neglect the states with t] = — 
which correspond to particles below the Dirac cone. The density of states at the tunneling 
energy eV is weighted by the probability density of the STM tip at position [x, y] for n=0. 
The contours for n = ±1 will be parametrized in terms of the polar angle 6 and transverse 
coordinate u. The proportionality factor J for the tunneling probability (not shown in the 
equation ) = JD(y; x,y) is a function of the distance between the tip and the sample. 
The notation D'^"'\V;x,y) represents the tunneling density for the different contours. 

IVB-The tunneling density of states D^^^^^V; x,y) for n = 

Summing up the single particle states weighted with occupation probability 
|t/i"=°'™'^'-\a:,2/)|2, we obtain a space dependent density of states for the two dimensional 
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boundary surface < x < ^ and the coordinate y is restricted to the regions | < y < f 
and ^ < y < |. We will perform the computation at the thermodynamic limit, namely 
we replace the discrete momentum j^k hj Y = and hj X = ^ where N = ^. We 
find for the dimensionless momentum q = qa the equations : q±{Y) = nV ± ;^tan~^[^|^] 
where 2tiX = pa = p. As a result we obtain the following density of states ^ 



dY TT q++P 

[ 



_lq-+f + j!P 



dY' TT + p2 

(50) 

Using this results, we compute the tunneling density of states in terms of the energy + eV 
measured with respect the chemical potential /i and the transverse energy e± = hvFq±- 

li(i + rT7^,A-(^)^ ) + i(i-rT7^,A-(:'" " " 



2' 7rL(// + T/)Y > + y' ' 2' nL{fjL + V)\j > + y' 
+^{H[fi + V-^]-H[fi + eV- E^,,]) ■ {{cos['J^-^—^y - -]f + (cos[^^^^— ^y - 

= (- — )^(-— )^ — , e - / dej ^ 



'hvp' L 4^G{x,y) '2 L ^ 2L 

for \y\ > d 



H[ij, + eV — ^jf-] is the step function which is one for fx + eV — ^jf- > and zero otherwise, 
a = ^ is the short distance cut-off and E^ax — hvpA < O.SeV is the maximal energy which 
restricts the vahdity of the Dirac model. We observe in the second line that the asymmetry 
in the density of states 1 ± \ lI^^^v) \J^ ~ (/Thf^ ) cancels. 

Equation (51) shows that the tunneling density of states is linear in the energy ji + eV (in 
the present case we have looked only for energies above the Dirac cone ). For the chemical 
potential fi = 120mV, the zero energy corresponds to the Voltage V = — 120my. The 



(51) 
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tunneling density of states has a constant part at energies ~ 0.2mV for —120mV < 
V < — 119.8mV. For V > — 119.8mV^ the density of states is proportional to /i + eV. 

In figure 2 we have plotted the tunneling density of states as a function of the coordinates 
X and y. The shape of the plot is governed by the the multiplicative factor e ^ '-x±iyi which 
governs the solutions in eq.(35). We observe that the density of state is maximal in the region 



\y 



< 105(2). 

Figure 3 shows the dependence on the voltage V and coordinate y. We observe the linear 
increase in the tunneling density of states which is maximal in the region \y\ < lO-B'-^''. 

IVC-The tunneling density of states D^^-^^^V, x,y;ri, .■r2M) for 2M dislocations. 

For many dislocations which satisfy Yl'^=i -B*-^''"-' = ( sum of the Burger vectors is zero 
) with the core centered at [x^,?/^] ,w = 1,2..2M the coordinate r = {x,y) — > [X{f),Y{r)] 
is replaced by [X{r) = x,Y{r) = y + ""^(frf^)]- Following the method used 

previously, we find the edge Hamiltonian with many dislocations takes the form: 

. 2M 

h"^^\w = 1,2...2M) ^ia'ldy - '-J^^^'^B^^'^^S^f - f^)] - la^d, (52) 

10 = 1 



jr{n=0,w=l,2...2M)y N TT --^(7 rr"? 7)l7{"=0)/ N 

'[x,y)c<. II e y{x-x^)+z(y-y^)> {x,y) 

.(n=0,i«=l,2...2M). TT -^( ,._..„,i,;,„_„.„ J ^(^=0) , 



As a result, the wave functions are given by: 

_ S(2) 

«;=1,2...2M 

^|n=0,i«=l,2...2M)^^^^^ ^ -Q )j.|n=0) ^53^ 

w=l,2...2M 

Using these wave functions, we find that the tunneling density of states is given by: 

10=1,2. ..2A/ 

In figure 4 we show the tunneling density of states for an even number of dislocations in the 
y directions which have the core on the ?/ = 0], w = 1,2,3,...2M). 

When the coordinate of the w = 1, 2, 3, ...2M dislocations is replaced by a continuum variable 
w which can be described by a domain a wall model: fi<io"^<^^'^-^<^'-'-(^x,y) = hvpl—icr^dy + 



, -j—f -B(2) , (x-xiu) . 

D^'^-"' {V,x,y;f I,.. f 2m) oc || e »)"+(y-j/™)"+a^^ (54) 



ia^d^ - a^KM{y)] where M{y) = sgn[y]\M{y)\ [25j. 

Using this model find that the tunneling density of states density Dd.omain-waiiiy.^^^y^ 
confined to |?/| < (the width W depends on the explicit form of the domain wall function 
M{y) and strength k) is given by: D'^°""'^''-'"''^\V ; x , y) oc (7^)^6-2'^ io"' <^y''''^^y'\ This show 
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the similarity between the result obtain from the domain— wall model and the large numbers 
of of dislocations given in equation (54). 

IVD-The tunneling density of states D'^'"'^^^\V;9,u) for the n = ±1 contours. 

Following the same procedure as used for the n — and using the eigenf unctions for 
n — ±1 we find : 

£,(-=±i)(y; e, u) = L>("=±^)(/i, V; 9, u)even + D^^=^'\f,, V; 9, u, f^Ud (55) 
For the even /c's, we solve for the momentum (5+ and and find: 



f R(2)'|2 oo / 7 _|_ 1 

27rRg{l)D{l)^yG{9,u) V ^aiV 

1 1 
2(-l)^ sinfQrii + 2'*(^' Qr)] sin[QrU - -k{1, Qr)] ■ 

(cos[/(^ + — ) - C(^, u){-r]{u) + cos[29])] - cos[(/ + 1)(^ + — ) + C(^, u){-r]{u) + cos[2^])]]; 

(56) 

Similarly for the odd /c's we find: 



27ri?5(l)i?(l)VG'(^,M)Q^J^Q_t^ V ^3^1) 

^^^-2a9,u)sH29] ^ e2^('^'")^H2^^1)((cOs[g,Ii - ^Ac(/, Qr)]f + {cOs[QrU + ^k{1, QrW) + 
2(-l)' COs[QrU + ^k(/, Qr)] COs[QrU - ^fi;(/, Qr)] ■ 

(cos[i(^ + — ) - C(^, u){-rj{u) + cos[29])] - cos[{l + 1){9 + —) + ({9, u){-rj{u) + cos[29])]] 



(57) 



For the present case the energy scale of the excitations is governed by the radius Rg{l) and 
width D. The spectrum is discrete and we can't replace it by a continuum density of states 
as we did for the case n = 0. 

In figure 5 we show the tunneling density of states at a fixed polar angle = | as a 
function of the voltage V. We observe that the density of states is dominated by high 
energy eigenvalues. This solutions are localized in energy. The range of the spectrum is 
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above ji + eV > 200my which is well separated from the low energy spectrum controlled by 
the n = contour (which ranges from —12{)mV to IQmV). 

Figure 6 shows the tunneling density of states as a function of the polar angle 6 for a 
fixed energy . The periodicity in 9 is controlled by the discrete energy eigenvalues. 

In figure 7 we show the tunneling density of states at a fixed voltage as a function of 
the polar angle < ^ < tt and width \u\ < 0.1. 

V-The charge current-the in plane spin on the surface of the /i^^ Hamiltonian 
A-The current in the absence of the edge dislocation for the h^^^^ 

Prom the Hamiltonian given in equation 1 we compute the equation of motion for the 
velocity operator: ^ — ^[x,h] — vpa^ , ^ — j^[y,h] — —vpcr^- We multiply the velocity 
operator by the charge (— e) and identify the charge current operators : Jx — {—e)vFcr^, 
Jy = {—e){—VF)a^. This also represent the '"real"' spin on the surface. Therefore, the 
charge current is a measure of the in-plane spin on the surface. 

Integrating over the y coordinate we obtain the current Jj'^' in the x direction. Using 
the eigenstates rp^q{x,y) and r^p^q{x,y) of the h^'^' Hamiltonian 



we find {rp^q{x,y)){a'^){rp^q{x,y)) = -(r_p,g(x, y))((7^)(r_p,g(x, y)) therefore, we conclude 
that the current /J-^- = is zero. 

VB-The current in the presence of the edge dislocation 

We will compute the current in the presence of the edge dislocation. The current operator 



Jf'^%x,y) = {-e)vF[a^el-a'el] = (-e) W - (-eK^( ^ , ) ^ {-e)vFa^ (58) 



quantization form for the current density. The operator is defined with respect the to 





J^'^^^{x,y) will be given in terms of the transformed currents. We find that the current 
density operator J^'^^'^{x,y) is given by: 
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shifted ground state |/i >= |0 > with the energy E — e — /i measured with respect the 
chemical potential and spinor field 'i>n=o{x,y)- 

Jf''{x,y) =< pi\<lfUo{x:y)Jf''{x,y)<lfn=o{^:y)\pi > (59) 

Using the spinor eigenfunction given in equation (35) and the second quantized form with 
the electron like operators aE,R,0(E,L and hole like (3e,r,(3e,l we find : 



The current is a sum of two terms computed with the eigen spinor obtained 
in equation (35): [C/{"=°'''^(x, y), C/f =°'''^(x, y)]^a2[C/{"=°'^)(x, y), ^^^^ and 
\U^^~^"^\x^y)^U^^~^"^\x,y)Ya'^\U^^~'^'^\x,y)^U^^~'^"^^^ which have opposite signs. 
Due to the parity violation caused by the dislocation, the density of states is asymmetric 
1 ± ~ L(^J^v) \J'^ ~ ^•]^tv'^'^ ) resulting in a finite current. We integrate over the transversal 
direction y and obtain the edge current /^=0'edge 



jn=,,ea,e ^ {-e)vE J_l ^ J] dy < fi\jf^^{x,y)\fi >- 



2 2 



-e)vF L .2/^ \ dx dy e ^^^tpt^> r r r — - — (hvp/L) ■ e\\ 



/'^ dx f '2 dy e " ^x'^+y'^+a'^ ' f f / 



2 2 



H[ijl — ^/(eii)^ + {^±y ] is the step function which is one for ^/(eii)^ + (e±)^ < A*- The 
single particle energies are e± — hvFq± and e|| — hvpp. For L 10~^m, chemical potential 
/X = 120mV and ~ 100 we find that the current /_j^'=o,e'^s'e -g -^-^ ^-j-^^ range of mA. 

To conclude, we have shown that the presence of an edge dislocation gives rise to a 
non-zero current which is a manifestation of the in-plane component of the spin on the two 
dimensional surface . Therefore a nonzero value 7^=0'e(ifle _^ q indication of the 

presence of the edge dislocation. This effect might be measured using a coated tip with 
magnetic material used by the technique of Magnetic Force Microscopy. 
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(61) 



Vl-Conclusions 

We have used the coordinate transformation method to investigate TI in the presence of 
deformations. We have computed the spin connection and the metric tensor for the three 
dimensional TI. This theory has been apphed to the surface of a TI with an edge dislocation. 
We have shown that the tunneling density of states is confined to two dimensional region 
n = and to high energy circular contours with n = ±1. The edge dislocations violate the 
parity symmetry. As a result a current which is a manifestation of in plane spin orientation 
is generated. The in plane spin orientation is a manifestation of the parity violation induced 
by the edge dislocation. We propose that scanning tunnehng methods might be able to 
verify our prediction. 
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Appendix -A 



We consider that a two dimensional manifold with a mapping from the curved space X", 
a = 1,2, to the local flat space x^, fi = x,y exists. We introduce the tangent vector [31| 
^/l(^) — ^^Jm^' = which satisfies the orthonormality relation e^{x)e^^{x) = 5a,b (here 
we use the convention that we sum over indices which appear twice). The metric tensor for 
the curved space is given in terms of the fiat metric 5a,b and the scalar product of the tangent 
vectors: e^(a;)ej!;(a;) = gn^i,{x). The linear connection is determined by the Christoffel tensor 



"^aA = -r;,.5A (62) 
The Christoffel tensor is constructed from the metric tensor gfj„u{x). 

^i-^ = ~\Y1 9^'^(^)l9'^9uA^) + df,guA^) - drQp^A^)] (63) 

T=x,y 

Next, we introduce the vector field V = V^da = V^d^ where a = 1, 2 are the components 
in the curved space and fi = x,y represents the coordinate in the fixed cartesian frame. The 
covariant derivative of the vector field V"' is determined by the spin connection o;^^ which 
needs to be computed: 



D^Vix) 



d,V''ix)+Lo^^,V'' 



(64) 



For a two component spinor, we can identify the spin connection in the following way: 
The spinor in the the curved space (generated by the dislocation) is represented by ^E'(X) 
and in the Cartesian space it is given by is given by \l'(x) 38|. The two component spinor 
represents a chiral fermion which transform under spatial rotation as spin half fermion: 



^(X) = e^'^^-2'"«*(f) 

(65) 
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We have used the anti symmetric property of the rotation matrix uja,b = —^b,a, and the 
representation of the generator S'^"'' in terms of the Pauh matrices. 
Therefore for a two component spinor we obtain the connection: 

D.^ix) = id, + ^a;^-''S,,,)vI/(f) ^ (d, + ^c^^'^ia,, a,])^ix) (66) 

Next we will compute the spin connection u'^'^ using the Christoffel tensor. In the phys- 
ical coordinate basis the covariant derivative D,V^{x) is determined by the Christoffel 
tensor: 

D,V^{x) = d,V'^{x) + TlX (67) 

The relation between the spin connection and the linear connection can be obtained from 
the fact that the two covariant derivative of the vector V are equivalent. 

D^V' = elD^V^ (68) 
Since we have the relation = e'^V" it follows from the last equation 

D,[el] = D^d^e'' = {D,d,)e^ + d,{D,e'') = (69) 

Using the definition of the Christoffel index and the differential geometry relation Vg^ di, = 
— r^^^c^A [31i], we obtain the relation between the spin connection and the linear connection: 

D,K] = d.elix) - ri,,e^(x) + co^^^eUx) ^ (70) 

Solving this equation, we obtain the spin connection given in terms of the Burger vector. 
We multiply from left equation (70) by the tangent vector and replace ^ with the 
representation given in equation (63). We use the metric tensor relations e^(x)e|^(x) = 6a^b, 
e^^{x)el{x) = g^,,yix). and find jsil: 

AeP'^e^'\d,e,,,-d^e,,,)el (71) 
We notice the asymmetry between e*^'" and Ca^v'- 
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For our case we have a two component the spin connection u;^^ and oUy"^ 
1 1 1 



(72) 

These equations are further simphfied with the help of equations (11 — 17) with = , 
el = 1 and the Burger tensor dxCy — dye^ = 



21^ 



l5(2)5(2)(f)U^^ - (1 - .^^)2..^^.., _ ^ 

l5(2)<5(2)(r-)[-^^] 
2 ^ 27r y2 + a;2J 

(73) 



and 



^^-'\i[a,e2 - a,e2] - ^^^'^e^[a,ei - a^ej] - \g^'^ e\g'^^^el{d,el - 9.e^]e^ = 

R(2) R(2) d(2) 

(74) 

To first order first the Burger vector 5 (2) the spin connections are given by : ux^ = —oj"^ ~ 
and ujf = -ujf ^ -\B^^)5\r). 
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y(x 



)/B 




FIG. 1: The contours {x{s)f + {y{s)-^)'^ = {^f for n = ±1, ±2, ±3(in decreasing size ),i?g(n) = 
n = corresponds to the equation y{s) = and |y| > d (see the text). The the distance is 
measured in units of the Burger vector B^'^\ 
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FIG. 2: The tunneling density of states forn = , ^ cx Z)("=°)(;g^, fi = UOmV). The right 
corner represents the intersection of the x coordinate which runs from 30 (right corner) to —30 and 
the y coordinate which runs from —30 (right corner) to 30 in units of the Burger vector. 
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D(x/B=-2;y/B,V;//=120) 



20 



FIG. 3: The tunneling density of states for n = as a function of y and V ^ D^"" ^\^2) — 
—2, -^(2)] ^ = 120mV). The voltage range is —120 < y < 50 and the y coordinate is in the range 
-30 < ^ < 30. 
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FIG. 4: Many Dislocations - with the core of the dislocations at [xw,y = 0] , w = 1,2...2M; 
The maximum of the tunneling density of states is confined along y = 0. The coordinates of the 
tunneling density of states are restricted to : —40 < < 40 and —20 < < 20. 
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D(6'=Pi/2,u/B=0.01,)U=120 ;V) 
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FIG. 5: The discrete tunneling density of states for n = 1, as a function of the voltage V 
Di-=^){V;e = ^,^,fi = 120mV) 
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D(6';u/B=0.01,V=280) 




FIG. 6: The tunneling density of states as a function of 6 D^'^=^\e] ^ = 0.01, V = 280mV, fi = 
120mV) 
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